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CROSS-REFERENCES TO RELATED APPLIC ATIONS 

This application claims priority from United States Provisional Paiem Application Ser. 
No. 60/406,998 filed on 29 August 2002. 

5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0001] The invention is in the field of nuclear magnetic resonance detection. More 
specifically, the invention is a method for distinguishing resonance signals within a 
mixture of substances with different transversal relaxation rates. 

10 

2. Background of the Art 

[0002] NMR relaxation measurements typically employ a static magnetic field to 
magnetize nuclei into an equilibrium state and an RF magnetic fieLi orthogonal to the 
static magnetic field to disturb this equilibrium state of nuclear magnetization. The RF 

15 magnetic field is typically applied in a form of short pulses that giv 3 rise to free induction 
decay signals (FED) in a nearby NMR antenna. There are two type* of relaxation 
processes related to energy flow and loss of coherence of spin magnetization in a 
population of nuclei. These relaxation processes are characterized )y two relaxation 
times: T p the spin-lattice relaxation time, and T v the spin-spin relaxation time. The 

20 spin-lattice relaxation time is the time constant associated with the ceturn of the 

longitudinal component to its thermal equilibrium state, which is piirallel to the static 
magnetic field. The spin-spin relaxation time is the time constant associated with the 
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decay of the transversal component of nuclear magnetization to zero. If che static 
magnetic field is in Z-direction of Cartesian coordinates, then the tiansversal component, 
inducing signal in the NMR antenna is in the X-Y plane. The net tiuclear magnetization 
in the X-Y plane decays to zero due to reversible and irreversible processes of de-phasing 
5 of spin isochromats. The reversible process is a result of macroscopic inhomogeneity of 

the static magnetic field. This process can be refocused into a spin echo signal by 
employing a refocusing RF pulse. 

[0003] A standard sequence of RF pulses used to measure the true iransversal relaxation 
10 (not related to the macroscopic inhomogeneity of the static magnetic field) is the CPMG 
sequence. The CPMG sequence is described, for instance, in Hxpe|imental Pulse NMR: 
A Nuts and Bolts Approach by Fukushima and Roeder. This sequence comprises a first 
excitation RF pulse that tilts the magnetization into X-Y plane followed by a plurality of 
refocusing RF pulses with the carrier frequency phase shifted by 90" with respect to the 
IS excitation pulse. The period of repetition of the plurality of the reft musing pulses is twice 

the length of time between the center of the excitation pulse and the center of the first 
refocusing pulse. The spin echo signal, which results from refocusing the spin 
isochromats, appears between refocusing pulses. The amplitudes o the echoes represent 
points on a T 2 relaxation curve. This curve is then decomposable into exponential teims 
20 in order to differentiate between the types of substances present ant/or, in the case of a 
fluid trapped in a porous stnicture, to characterize the porous medi* , In a mixture of 
substances, though, it is almost impossible to differentiate between substances just by the 
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exponential decomposition if the substance with the shorter decay rate is of very small 
quantity compared to the others. A systematic error in the relaxati :>n curve measurement 
caused by non-ideality of the CPMG pulse sequence and inhomog meity of the static 
magnetic field will most likely exceed the signal from the short relaxation substance. 

5 

[0004] Another standard NMR technique for measuring parameters of different 
substances is based on analyzing the FID. In the case of a high ho nogeneity of the static 
magnetic field (typically a large and expensive magnet is required and a short dead time 
of the NMR system, a high precision of separation between substa ices with short and 
10 long transversal relaxation can be achieved. An appropriate analysis of the FID to 

improve accuracy of the separation is described, for example, in U.S. Patent 5,519,319 
issued to Smith et al. In many practical situations, the FID signal i ; defined by the 
inhomogeneiry of the static magnetic field. No differentiation between substances is 
possible based on FID in this case. 

15 

[0005] Another known method for measurement of a substance with short relaxation in a 
mixture is described in U.S. Patent 6,232,778 issued to Speier et a . The method is based 
on repeating the NMR experiment quickly enough in order to prevsm the magnetization 
vector of a long NMR relaxation substance to recover its longitudi nal component and 
20 therefore effectively preventing participation in the accumulated si gnal. Following this 
procedure, a multi-exponential data analysis is performed in order to filter out the 
residuals of the long relaxation component. 
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[0006] Another way to increase the precision of NMR based-diffeientiation beiween two 
substances in a mixture is described in U.S. Patent 6,091,242 issued to Hanawa. A first 
RF pulse is applied to the sample in order to invert the equilibrium magnetization. For 
both substances, the longitudinal component of the nuclear magneiization starts to evolve 
5 toward the equilibrium. Due to difference in the longitudinal rela?i ation times, it is 

possible to choose an appropriate wait time after the application ol the inversion pulse to 
get a zero crossing on the longitudinal magneiization curve for om: of the substances in 
the mixture. An excitation pulse applied at the moment of the zen crossing of one 
substance will not cause any signal for that substance. As a result, even small amount of 
10 the other substance can be effectively detected. This method iequi res prior knowledge of 
T\ or T| distribution for at least one of the substances. 

[0007] The methods of Speier 778 and Hanawa '242 are ineffecti - 'e if the substances 
with different transversal relaxation times have comparable longitudinal relaxation times, 

15 This situation occurs in many practical cases. One example is relaied to detection of a 
small amount of a heavy oil, for example, Herro Negro, having N\tR relaxation 
parameters T 2 = 0.5ms and T t = 40 ms, in the presence of lighter oii Is or water in a porous 
media. The substances with longer relaxation rates may have the relaxation limes T { = T 2 
in the range 20-100 ms that are comparable with the T x of the heav / oil. Another 

20 example is found in detecting bone bound water in the human body. The bound water T 2 
relaxation time in the frequency range 2-10 MHz is typically 0.3-0. 7 ms while 
longitudinal relaxation time T, is 20-30 ms. The other fluids in the body have T 2 and T| 
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NMR relaxation times in the range 10-200 ms. 

[0008] Neither shortening the rime between repetitions nor selective T, relaxation can be 
effective in these cases. Both known techniques will cause the san e extent of NMR 
5 signal reduction for long T 2 relaxation components as for the short T 2 relaxation 

component, 

[0009] There is a need for a new pulse sequence that enables elimination of the signal 
from a long relaxation substance without spoiling the relaxation signal from a short 
10 relaxation substance, even if the two substances have the same lonptudinal relaxation 
times. The invention disclosed herein addresses this need. 

SUMMARY OF THE INVENTION 

15 [0010] The present invention is an improved method of measuring i small quantity of a 
first substance, said substance having a shojr T 2 relaxation time, in the presence of a 
second substance, said second substance having a long T 2 relaxation time. In a static 
magnetic field, a first excitation RF pulse is applied to the mixture substances. An 
inversion pulse is applied to refocus the nuclear magnetization of the long relaxation 

20 substance. At the moment of maximum spin echo of the long relax ation substance, a 
second excitation pulse is applied. At the moment immediately aftur the application of 
the second excitation pulse, the longitudinal component of the long relaxation substance 
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is inverted from its equilibrium, while che longitudinal component of the short relaxation 
substance is substantially zero. 

[0011] By observing the longitudinal component of the long relaxei tion substance alone, it 
5 is possible to obtain a time after being inverted by the second exciution pulse value at 
. which said longitudinal component of said long relaxation substance will cross zero. At 
the moment of zero crossing, a third excitation pulses is applied fo lowed by a sequence 
of inversion pulses. The combination of said third excitation pulses and the sequence of 
inversion pulses comprise a CPMG pulse sequence. If it is not posii ible to determine zero 

10 crossing by the method stated above, it can be noted that applying 1 he CPMG sequence 

produces a set of transverse spin echoes for non-zero longitudinal t omponents of the long 
relaxation substance. Zero crossing can thus be determined by min imization of said spin 
echoes. When said third excitation pulse is applied at the moment of said zero crossing, 
the signal of the long relaxation substance is minimized and the rel ixation signal of the 

15 shon relaxation substance can then be observed. 

BRIEF DESCRIPTION OF THE DRAWIN GS 

[0012] The present invention is best understood by reference to the; accompanying figures 
20 in which: 

FIG. 1 (Prior Art) depicts a typical NMR logging tool conveyed within a borehole. 
FIG. 1A (prior art) shows the magnet and antenna configuration o " a preferred device for 
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the method of the present invention. 

FIG- 2A shows the pulse sequence employed in the present invent on. 
FIG. 2B shows the evolution of nuclear magnetization component*; for a long T, and TV 
relaxation substance to the pulse sequence of Fig. 1. 
5 FIG. 3 shows the evolution of nuclear magnetization components lor a short T 2 
relaxation substance to the pulse sequence of Fig. 1. 

FIG- 4A shows a computer simulation of the evolution of nuclear Magnetization 
components in a mixture of short T 2 and long TVrelaxation substan ;es> 
FIG. 4B shows an expanded view of Fig. 4A. 
10 FIG. 4C shows an expanded view of Fig. 4A with reduced waictin:e. 

FIG. 4D shows an expanded view of Fig. 4A with expanded wait ti me, 

FIG. 5A shows the result of modeling corresponding to a short T 2 relaxation substance. 

FIG- SB shows an expanded view of Fig. 5A. 

15 DESCRIPTION OF THE PREFERRED EMBOD IMENT 

[0013] Figure 1 depicts a borehole 10 which has been drilled in a typical fashion into a 
subsurface geological formation 12 to be investigated for potential hydrocarbon- 
producing reservoirs. An NMR logging tool 14 has been lowered u,to the hole 10 by 
means of a cable 16 and appropriate surface equipment represented diagrammatically by a 

20 reel 18 and is being raised through the formation 12 comprising a pi urality of layers 12a 

through 12g of differing composition, to log one or more of the fon nation's 
characteristics. The NMR logging tool 23 is provided with bowspri ngs 22 to maintain the 
584-28566US 8 
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tool in an eccentric position within the borehole with one side of the tool in proximity to 
• the borehole wall. Signals generated by the tool 23are passed to the surface through the 
cable 16 and from the cable 16 through another line 19 to appropriate surface equipment 
20 for processing, recording and/or display or for transmission to a lother site for 
5 processing, recording and/or display. It is to be noted that processi ig may also be done, 
in part or in whole, by a downhole processor. 

[0014] Fig- 1A schematically illustrates a preferred embodiment o cthe present invention 
for the structure of ihe NMR tool This is described fully in US Parent 6348792 to Beard 

10 et al, having the same assignee as the present invention and the con tents of which are 
fully incorporated herein by reference. It is to be noted that the me hod of the present 
invention would be equally effective with other prior ait devices usi ed for NMR 
measurements in boreholes. The tool cross-sectional view in Fig. IA illustrates a main 
magnet 67, a second magnet 68, and a transceiver antenna, compri ring wires 69 and core 

15 material 60. The arrows 71 and 73 depict the polarization (e,g., frc tn the South pole to 
the North pole) of the main magnet 67 and the secondary magnet ( 8. A noteworthy 
feature of the arrangement shown in Fig, 1A is that the polarization of the magnets 
providing the static field is towards the side of the tool, rather than towards the front of 
the tool (the right side of Fig. 1 A.) as in prior art devices. The imj lortance of this rotated 

20 configuration is discussed below. 

[0015] The second magnet 68 is positioned to augment the shape o f the static magnetic 
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field by adding a second magnetic dipole in close proximity to the IF dipole defined by 
the wires 69 and the soft magnetic core 60. This moves the center of the effective static 
dipole closer to the RF dipole, thereby increasing the azimuthal exisnt of the region of 
examination* the desirability of which has been discussed above. The second magnet 68 
5 also reduces the shunting effect of the high permeability magnetic «: ore 60 on the main 
magnet 67: in the absence of the second magnet, the DC field woul i be effectively 
shotted by the core 60. Thus, the second magnet, besides acting as a shaping magnet for 
shaping the static field to the front of the tool (the side of the main nagnet) also acts as a 
bucking magnet with respect to the static field in the core 60. Thosi: versed in the art 
10 would recognize that the bucking function and a limited shaping co Jld be accomplished 
simply by having a gap in the core; however, since some kind of field shaping is required 
on the front side of the tool, in a preferred embodiment of the invention, the second 
magnet serves both for field shaping and for bucking. If the static f eld in the core 60 is 
close to zero, then the magnetosrrictive ringing from ihe core is sub itantially eliminated 

15 

[0016] A novel RF pulse sequence for performing NMR relaxation measurements is 
presented in Figure 2A. The sequence is written as 
A 1 -t-B 1 -t-A2-TW^ A3^B2i ( 1 ) 

The sequence comprises a first excitation pulse Al (101), a first refi>custng pulse Bl 
20 (102) delayed by a time x from the first excitation pulse Al (101), a forced inversion 
pulse A2 (103) delayed by the time t from the first refocusing pulse B 1 (102), and a 
second excitation pulse A3 (104) delayed by the time TW from the breed inversion pulse 
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A2 (103). The sequence also comprises a plurality of refocusing j utses B2 t (105) 
subsequent to the second excitation pulse A3 (104). The second excitation pulse A3 
(104) and the plurality of refocusing pulses B2, (5) together form 21 CPMG pulse 
sequence. The length of time TW depends on the time constant of the longitudinal 
5 relaxation of the spins of protons contained in a substance having ; l long transverse 

relaxation rime. The excitation pulses Al (101) and A3 (104) and the forced inversion 
pulse A2 (103) are 90* pulses which tip the nuclear magnetization vector generally 
through an angle of 90*. The refocusing pulses Bl (102) and B2j ( 105) are 180" pulses, 
which tip the nuclear magnetization vector generally through an an gle of 180°. However, 
10 the B pulses could also have a reduced tipping angle (less than 18C °) as discussed in US 
Patent 6,163,153 to Reiderman et al. 

[0017] Figure 2b is an example of a computer simulation showing how the transversal 
and longitudinal nuclear magnetization components evolve for a si bstance with long 

15 relaxation times during application of the pulse sequence presentee in Equation 1. To 

illustrate the invention, a long relaxation substance is stated as hav ng T, * T 2 ' * 30msec 
where T 2 ' indicates a substance having a long relaxation time. A si atic magnetic field 
oriented in the +z-direction is originally applied to the substance, thereby aligning the 
individual nuclear magnetization vectors along the +z axis- The m ignetization vector of 

20 the substance, comprising the individual nuclear magnetization ve<: tors, is originally in an 
equilibrium state and is oriented in the +z direction. The first excii: ation pulse Al (201) 
rotates the nuclear magnetization vector from its equilibrium state (parallel to Z-axis) by 
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90" and into an X-Y plane having the z-axis as its normal. The subsequent magnetization 
vector component, stated as M x (210) without loss of generality, ch :n de-phases 
(corresponding to a FID) due mainly to the inhomogeneity of the s I atic magnetic field. 
First refocusing pulse Bl (202) rotates the dephased vectors by 18C ° enabling a spin echo 
5 to occur at the time 2t. Due to the long T, times, the longitudinal i; omponent (210) 

experiences insignificant growth and is substantially zero at time t=i2r. Refocused at time 
t=2x, the transversal magnetization M, (210) is then rotated 90° by the forced inversion 
pulse A2 (203) to lie in the direction opposite to the equilibrium, in the -z- direction. 
After the application of pulse A2 (203), the longitudinal magnetize ion (212) begins to 

10 evolve from its inverted state back to the equilibrium state. On its i ipproach towards 

equilibrium, (212) crosses through a point of zero net magnetization, herein referred 
co as the zero axis. At the moment of generally crossing zero, the second excitation pulse 
A3 (204) is applied. The longitudinal component M,, prior to applk arion of pulse A2 
(203) is tilted into the X-Y plane and becomes the transverse compi mem (214). This 

15 transverse component (214) is insignificant and diminishing for the time TW between 
pulse A2 (203) and pulse A3 (204). For this example, the applied pulse (204) results in 
an absence of any transverse component of the magnetization and, <; orrespondingly, no 
signal induced in the NMR antenna during the application of the plurality of refocusing 
pulses B2 { (205). 

20 

[0018] An essentially different situation takes place in the case of a substance having a 
short transverse relaxation rate, indicated as T 2 with Ti»T 2 , as is piesented in Figure 3. 
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By setting T according to the inequality T 2 « T <<Z X we ensure that no spin echo from the 
short T 2 relaxation substance is created^ At the moment immediate ly prior to application 
of the B 1 pulse (302) the transverse component has significantly d :phased to zero. As a 
consequence, pulse Bl does not produce a spin echo at time t ss 2t and immediately prior 
5 to the first inversion pulse (303), the transverse component is subs antially zero. Thus, 
after application of the first inversion pulse (303), the longitudinal component of the 
magnetization M, (312) is generally zero. This zero value is in coi tiast to the substantial 
inverted component of the long relaxation substance shown in Figure 2b. After the 
elapsed time TW and immediately prior to the application of the s<: cond excitation pulse 

10 (304), the magnetization vector M, (312) has practically reached it: , equilibrium. Again, 
this is in contrast to the situation with the long T, 1 relaxation subst ince of Figure 2b, 
wherein the longitudinal magnetization vector (212) is generall 1 zero at the moment 
of application of the pulse A3 (204). The pulse (304) hence tips thsMj (312) 
magnetization into the X-Y plane to generate a FID signal (314) in the NMR receiving 

IS antenna. 

[0019] Figures 4A, 4B present the result of computer modeling for a mixture of two 
substances having different relaxation rates and reacting to the infl ience of the pulse 
sequence of Equation 1. The slow relaxation substance has relaxation times T x = T 2 ' = 30 
20 msec. The fast relaxation substance, in quantity of 5% of the long relaxation substance, 
has relaxation parameters T, = 30 msec, T 2 = 0.5 msec. It is clear that for a mixture of 
long relaxation signals there will be a point after the application of the inversion pulse 
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(403A) at which the longitudinal component of magnetization (413 A) crosses the zero 
axis. Figure 4B presents an enlarged section of Figure 4A in the vicinity of the second 
excitation pulse A3 (404A) ? with pulse 404A of Figure 4B representing pulse 404A of 
Figure 4A. At the moment of zero crossing, che relatively small sijj jial of a short 
5 relaxation FID signal (414B) can be measured. The longitudinal ci: .mponent M 2 (416B) 
experiences growth toward equilibrium as well as 180° inversions ilue to the pulse 
sequence B2, (405B). 



[0020] The moment of application of pulse A3 (404A ) is determin id by observing the 
10 point at which the longitudinal component M, (412B) only for the long relaxation 

substance crosses the zero axis. In case the zero-crossing point cam not be determined in 
advance by using measurement results on the long relaxation subsfci Jices only, or if the 
moment of crossing the zero axis varies due to temperature dependence of T! relaxation 
time for the long relaxation substances, then the procedure described below can be 
15 preferably used to find the zero crossing point. The CPMG sequem :e, comprising pulse 
(404B) and the plurality of pulses (405B), generates a train of spin echoes. If the delay 
after pulses in the plurality of pulses (405B) exceeds the T 2 relaxati m time for the short 
relaxation substance, then any echoes seen between pulses (405B) represent only those 
echoes from the long relaxation substances. Thus, by varying the ti me interval TW 
20 between the pulses A2 (403A) and A3 (404A) in Figure 4A, it is possible to zero out the 
echo amplitudes M x (414B) and thus to use the zeroing of echo amplitudes as an 
indicator. This method enables an operator to obtain a high precisi: >n elimination of the 
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signal from the slowly relaxing substance. Figure 4B is an illustraiion for the situation of 
substantially zero echo amplitudes and, correspondingly, for the fa. t that the FID signal 
(414B) after the pulse (404B) represents the contributions significantly from the short 
relaxation substance. 

5 

[0021] To further illustrate this technique, Figures 4C and 4D represent the situations 
when TW is 8% smaller and 8% larger, respectively, than its value used in Figure 4B. 
Comparison of Figures 4C and 4D illustrates the high sensitivity of the echo signal M, to 
relatively small departures of TW from its optimum value. In Figti re 4C, the 

10 component (412C) is has a net negative value at the moment of ap] ilication of the second 
excitation pulse (404C). The resultant M, component (418C) subsequent leads to M, 
echoes (420C) having substantially negative polarity. Similarly, In Figure 4D, the IvI^ 
component (412D) has a net position value at the moment of application of the second 
excitation pulse (404D). The resultant M, component (418D) subsequent leads to M, 

1 5 echoes (420D) having substantially positive polarity. 

[0022] In another embodiment of the invention, after the wait time TW needed to zero out 
the free induction decay (FED) signal M, (414B) is determined, the TW is altered slightly. 
This is followed by a single A3 pulse. Analysis of the free inductic m decay signal can 
20 then be used for determining the relative quantities of the fast and • lowly relaxing 
components. 
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[0023] Figure 5A and SB shows the result of modeling correspond: ng to the 5% quantity 
of short T 2 relaxation substance only affected by the pulse sequeno: of Equation 1. 
Figure 5A shows the growth of the longitudinal component M,, (51 1) and the FID (520) 
that results from applied pulse 504. Figure 5B is an enlargement o : Figure 5A in the 
5 vicinity of pulse 504, relabeled as 504* in Figure 5B. Comparison >f the short relaxation 
FID (520') signal after the second excitation pulse (504*) in Figure 5B with the recovered 
FID signal 414 in Figure 4B shows that the signal from a large amc unt of the long 
relaxation substance in the mixture is substantially eliminated from Figure 4B. Thus, the 
proposed technique enables measuring a very small amount of the i\ hort relaxation 
10 substance mixed with a large amount of long relaxation substances 

[0024] The present invention has been described with reference to ;k wireline device. 
However, the principles of the invention may also be embodied in ::nd used with MWD 
devices conveyed on a drilling tubular such as a drillstring or coilec rubing. 



15 



[0023] While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the ait. It is intended 
that all variations within the scope and spirit of the appended claim ; be embraced by the 
foregoing disclosure. 
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